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This  report  details  the  application  of  a  two-dimensional,  laterally 
averaged,  water  quality  code,  CE-QUAL-W2,  to  the  Savannah  River  Estuary,  where 
resolution  of  both  vertical  and  longitudinal  gradients  was  desired.  The 
report  is  Intended  for  an  audience  that  is  familiar  with  numerical  water  qual¬ 
ity  modeling. 
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PART  I.  INTRODUCTION 


General  Comments 


1.  Coastal  construction  and  operation  activities  at  US  Army  Corps  of 
Engineer  (CE)  water  resource  projects  can  impact  the  quality  of  surface  estu¬ 
arine  waters.  Numerical  water  quality  models  are  required  to  assess  these 
impacts,  to  evaluate  various  structural  and  operational  alternatives  for  water 
quality  control,  and  to  determine  cause  and  effect  relationships  of  water 
quality  problems.  Frequently,  the  impacts  are  related  to  vertical  stratifica¬ 
tion,  such  as  bottom  anoxia  or  transport  of  larvae  and  contaminants.  Con¬ 
sideration  of  these  impacts  requires  the  use  of  numerical  water  quality  models 
with  resolution  in  both  the  vertical  and  longitudinal  dimension. 

2.  This  report  describes  the  application  of  CE-QUAL-W2  to  an  estuarine 
system  where  vertical  and  longitudinal  resolution  was  required.  CE-QUAL-W2  is 
a  two-dimensional,  laterally  averaged,  water  quality  model  that  was  developed 
to  assist  in  the  analysis  of  water  quality  problems  in  estuaries  and  reser¬ 
voirs  where  buoyancy  is  important  and  where  lateral  homogeneity  can  be 
assumed.  Detailed  description  and  user  guidance  are  provided  in  the  report 
"CE-QUAL-W2 :  A  Numerical  Two-Dimensional,  Laterally  Averaged  Model  of  Hydro¬ 
dynamics  and  Water  Quality;  User’s  Manual"  (Environmental  and  Hydraulics  Labo¬ 
ratories  1986). 

3.  This  application  consisted  of  two  levels  of  constituent  kinetic 
formulations:  (a)  dissolved  oxygen  (DO)  and  biochemical  oxygen  demand  (BOD) 
and  (b)  DO,  BOD,  nutrients,  and  algae.  Sensitivity  analyses  were  conducted 
using  both  levels  of  constituent  kinetic  formulations.  Tidal  boundary  DO  and 
BOD  concentrations  were  varied  in  the  DO-BOD  formulation,  and  the  system 
tldally  averaged  DO  was  noted  as  the  response  variable.  Tidal  boundary  algal 
concentration  and  maximum  algal  growth  rate  were  varied  in  the  second  formula¬ 
tion  with  system  tidally  averaged  algal  concentration  noted  as  the  response 
variable.  The  purpose  of  the  sensitivity  experiments  was  to  investigate  the 
Importance  of  tidal  boundary  constituent  concentrations.  Experiments  were 
conducted  using  the  DO-BOD  formulation  to  demonstrate  the  utility  of 
CE-QUAL-W2  for  investigating  changes  in  upstream  regulation  and  BOD  loadings. 
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4.  The  Savannah  River  basin  is  located  in  southwestern  North  Carolina, 
western  South  Carolina,  and  northern  Georgia.  The  confluence  of  the  Tugaloo 
and  Seneca  rivers  at  Hartwell,  Ga.,  form  the  Savannah  River,  which  flows  along 
the  Georgia-South  Carolina  state  line  to  the  Atlantic  Ocean.  Stream  dis¬ 
charges  in  the  basin  are  typically  high  in  the  winter  and  early  spring, 
decrease  during  the  summer,  and  remain  low  during  the  autumn.  Freshwater  flow 
is  regulated  through  upstream  reservoirs. 

5.  The  Savannah  River  is  tidal  for  approximately  70  km,  from  its  mouth 
to  near  Ebenezer  Landing  (River  Mile  (RM)  44.7).  A  deep-draft  harbor  is  main¬ 
tained  at  Savannah,  Ga.  The  channel  of  the  Savannah  Harbor  is  maintained  from 
the  Atlantic  Ocean  upestuary  to  RM  21.3.  The  tidal  variation  of  the  Savannah 
River  Estuary  of  over  2  m  is  exceeded  along  the  east  coast  only  by  the  estua¬ 
rine  waters  of  Maine. 

6.  A  unique  feature  of  the  harbor  designed  to  exploit  the  large  tidal 
variation  is  a  sediment  basin  and  tide  gate  to  reduce  sedimentation  within  the 
Savannah  Harbor  (Figure  1).  The  sediment  basin  tide  gate  plan  consists  of  a 
dredged  sediment  basin  at  the  lower  end  of  Back  River,  construction  of  a  tide 
gate  in  Back  River,  and  various  dredging  and  channel  Improvements.  The  tide 
gate  is  controlled  to  open  during  the  flood  portion  of  the  tidal  cycle  and 
then  to  close  at  high-water  slack.  The  tidal  prism  in  the  Back  and  Middle 
rivers  augment  ebb  flow  through  the  Front  River.  The  higher  resulting  ebb 
velocities  in  Front  River  tend  to  move  sediment  past  the  Savannah  Harbor  and 
concentrate  much  of  the  shoaling  in  the  sediment  basin. 


Objectives 

7.  The  primary  objective  of  the  Environmental  Impact  Research  Pro¬ 
gram  (EIRP)  work  unit  entitled  "Environmental  Impacts  of  Modifying  Estuarine 
Circulation  and  Transport  Processes"  was  the  development  of  a  CE  estuarine 
water  quality  modeling  capability  that  could  relate  changes  in  water  quality 
to  physical  changes,  such  as  circulation.  Capability  development  included  the 
selection  and  refinement  of  numerical  water  quality  codes.  The  numerical 
water  quality  model  CE-QUAL-W2  was  selected  for  estuaries  where  vertical  and 


VICINITY  MAP 

SCALE 

100  0  100  200  KM 


SEDIMENT  BASIN  -  TIDE  DATE 
STRUCTURE 


Figure  1.  Study  area  and  sediment  basin-tide  gate 

structure 


longitudinal  resolution  was  required  and  where  lateral  homogeneity  could  be 


assumed . 


8.  The  Savannah  River  Estuary  was  selected  for  the  test  application 
because  of  the  availability  of  data  and  information  from  previous  modeling 
studies.  The  US  Army  Engineer  Waterways  Experiment  Station  (WES)  coordinated 
a  comprehensive  water  quality  survey  of  the  Savannah  River  Estuary  during  the 


autumn  of  1979  and  late  summer  of  1980  (Pennington  and  Bond  1981;  Shingler  et 
al.  1981).  This  survey  was  conducted  to  assist  the  US  Army  Engineer  District, 
Savannah  (SAS) ,  In  the  preparation  of  a  plan  of  study  that  resulted  In  the 
application  of  the  one-dlmenslonal  water  quality  numerical  code  by  Lawler, 
Matusky,  and  Skelly  Engineers  (LMS)  (1983)  to  Investigate  the  DO  resources  of 
the  Savannah  River  Estuary.  In  addition  to  the  availability  of  data,  the 
Savannah  River  Estuary  Is  deep  and  narrow  so  that  the  lateral  homogeneity 
assumption  was  valid. 

9.  Additional  objectives  of  the  present  study  included  the  identifica¬ 
tion  of  minimum  data  requirements  and  recommended  sampling  designs  for  estu¬ 
arine  application  and  the  Identification  and  formulation  of  procedures  for 
data  summarization  and  preparation  for  model  application. 

10.  This  report  provides  the  results  of  the  application  of  CE-QUAL-W2 
to  the  Savannah  River  Estuary  for  the  period  August  1980.  Part  II  of  this 
report  describes  the  application  procedure,  Part  III  presents  the  study 
results,  and  Part  IV  presents  the  study  conclusions. 

Water  Quality  Data  Resources 

11.  The  survey  results  (Pennington  and  Bond  1981;  Shingler  et  al.  1981) 
and  the  effluent  data  reported  In  LMS  (1983)  were  used  In  this  study.  Table  1 
summarizes  the  survey  data  used  for  the  Savannah  River  application,  as 
follows: 

a.  The  center-line  survey  was  conducted  weekly  during  August  1980 
at  16  stations  along  the  main  channel  of  the  estuary  and 
provided  water  quality  constituent  profiles  used  for  model 
adjustment. 

b.  The  13-hr  survey  was  conducted  four  times  at  weekly  intervals 
at  the  downstream  study  boundary.  Each  weekly  study  consisted 
of  water  quality  constituent  profiles  collected  hourly  over  a 
13-hr  period.  The  results  of  the  13-hr  survey  were  used  to 
specify  tidal  conditions. 

c.  The  dlel  study  was  conducted  weekly  at  the  upstream  boundary  of 
the  study  area.  Surface  to  bottom  profiles  were  taken  every 

4  hr  for  measurements  of  temperature,  DO,  and  conductivity. 

The  results  of  the  dlel  study  were  used  to  specify  freshwater 
flow  boundary  conditions. 


Survey  Data  Used  for  the  Savannah  Rlv«r  Estuary  Application 
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d.  Experiments  to  estimate  production  by  the  light-dark  bottle 
method  were  conducted  weekly  at  the  downstream  boundary  con¬ 
comitant  with  the  13-hr  survey.  The  results  of  the  productiv¬ 
ity  survey  were  used  to  estimate  the  light  attenuation 
coefficient  and  the  fraction  of  incident  solar  radiation 
absorbed  at  the  water  surface. 

e.  Sediment  oxygen  demand  rates  were  determined  at  three  locations 
along  the  lower  reach  of  the  study  area.  The  results  were 
extrapolated  to  assign  sediment  oxygen  demand  rates  to  the 
entire  study  area. 

J[.  The  30-day  hydrodynamic  survey  provided  water  surface  level 
records  that  were  used  to  specify  the  tidal  boundary  elevation 
and  were  also  used  for  model  adjustment. 

g.  The  13-hr  hydrodynamic  study  consisted  of  tidal  velocity  mea¬ 
surements  taken  at  1-hr  intervals  at  the  upstream,  downstream, 
and  two  intermediate  locations  on  14  August  1980.  The  results 
were  used  for  model  adjustment. 


PART  II.  APPLICATION  PROCEDURE 


Geometric  Schematlzatlon 

12.  The  modeled  study  area  extended  from  the  confluence  of  the  Atlantic 
Intracoaatal  Waterway  (RM  5.9)  to  Ebenezer  Landing  (RM  44.7)  and  included  the 
reentrant  Back  River  connected  with  the  main  branch  at  RM  10.6  and  RM  19.8 
(Figure  1). 

13.  The  Savannah  River  Estuary  Is  represented  In  the  model  as  three 
branches:  Branch  1  was  the  main  channel  between  RM  5.9  and  RM  44.7;  Branch  2 
was  the  upstream  section  of  Back  River  from  the  tide  gate  to  its  intersection 
with  the  main  channel  at  RM  19.8;  and  Branch  3  was  the  downstream  section  of 
Back  River  from  the  tide  gate  to  Its  intersection  with  the  main  channel  at 

RM  10.6.  The  tide  gate  operation  was  represented  as  withdrawals  from  Branch  3 
during  flood  tides  with  Insertion  of  the  withdrawals  as  flows  Into  Branch  2. 

14.  Branch  1  Ax  was  2  km,  and,  the  Ax  for  Branches  2  and  3  was  1  km, 
where  Ax  Is  the  cell  length.  Branch  1  Ax  of  2  km  was  a  balance  between  suf¬ 
ficient  refinement  to  characterize  phenomena  of  interest  (effluent  discharges, 
branch  Intersections,  and  prototype  data  sampling  stations)  and  unnecessary 
computational  costs.  The  selection  of  Ax  of  1  km  for  Branches  2  and  3  was 
mandated  by  the  requirement  for  a  minimum  of  three  computational  segments  per 
layer.  The  plan  view  segmentation  with  segment  numbering  Is  displayed  In 
Figure  2.  The  main  channel  consisted  of  Segments  2  through  33;  the  upstream 
section  of  Back  River  consisted  of  Segments  36  through  43;  and  the  downstream 
section  of  Back  River  consisted  of  Segments  46  through  49. 

15.  One  study  objective  was  to  characterize  the  vertical  distribution 

3 

of  water  quality  constituents.  At  173-m  /sec  discharge  of  the  Savannah  River 
at  New  Savannah  Bluff  Dam,  the  surface  elevation  at  Ebenezer  Landing  is 
3.47  m.  The  surface  elevation  at  the  tidal  boundary  (RM  5.9)  at  mean  low 
water  (MLW)  Is  -1.04  m.  The  resultant  surface  slope  Is  7.2  *  10  ^  resulting 
In  a  total  elevation  difference  of  4.5  m  between  the  model  boundaries.  Initial 
simulations  Indicated  that  Instantaneous  water  surface  elevations  extending 
over  several  model  layers  resulted  in  numerical  Instabilities.  Therefore,  a 
surface  layer  thickness  of  5  m  was  selected,  with  the  remaining  layers  1  m 
thick.  The  selection  of  varying  layer  thicknesses  required  modification  of 


Figure  2.  Plan  view  segmentation  and  segment 
numbering 
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the  code.  Figure  3  displays  the  vertical  and  longitudinal  geometric 
schematlzatlon . 

16.  The  procedure  used  to  obtain  cell  widths  involved  the  use  of  the 
Hydrologic  Engineering  Center  code  GEDA  and  the  preprocessor  GIN.  GEDA  uses 
cross-sectional  data  from  ranges  at  irregular  Intervals  and  produces  estuary 
widths  as  functions  of  elevation  at  regular  Intervals.  GIN  uses  GEDA  output 
as  input  and  produces  the  geometry  (BA)  cards  read  by  CE-QUAL-W2. 

17.  Estuary  geometry  data  were  abstracted  from  the  following  sources: 

a.  Natural  Chart  11512,  Savannah  River  and  Uassaw  Sound,  National 
Oceanic  and  Atmospheric  Administration  (NOAA) ,  16  July  1977. 

b.  Nautical  Chart  11514,  Savannah  River,  Savannah  to  Brier  Creek, 
South  Carolina-Georgia,  NOAA,  7  January  1978. 

c.  Annual  Survey-1981,  Savannah  Harbor,  Georgia,  US  Coastal  High¬ 
way  No.  17  to  Sea,  US  Army  Engineer  District,  Savannah. 

d.  Navigation  Charts,  Savannah  River,  Georgia  and  South  Carolina, 
Savannah  to  Augusta,  June  1980,  US  Army  Engineer  District, 
Savannah. 

e.  Design  Memorandum  2,  Sediment  Basin,  Savannah  Harbor,  Georgia, 
Tide  Gate  Structure,  Vol  I  of  II,  US  Army  Engineer  District, 
Savannah,  31  January  1969. 

f .  Design  Memorandum  3,  Sediment  Basin,  Savannah  Harbor,  Georgia, 
Fresh  Water  Control  Plan  and  Appurtenances,  US  Army  Engineer 
District,  Savannah,  3  April  1970. 

18.  The  geometry  source  publications  related  depths  in  feet  to  MLW  or  a 

low-water  plane  above  Big  Collins  Creek  corresponding  to  a  discharge  of 
3 

6,100  cfs  (173  m  /sec)  at  New  Savannah  Bluff  Dam.  The  following  values  were 
used  to  relate  depth  to  0.0  ft*  MLW: 


RM 

MSL,**  ft 

5.9 

-3.4 

7.7 

-3.5 

9.4 

-3.6 

10.9 

-3.6 

12.8 

-3.6 

14.4 

-3.6 

16.3 

-3.6 

18.0 

-3.6 

19.6 

-3.6 

21.5 

-3.5 

*  To  convert  feet  into  metres,  multiply  by  0.3048. 

**  MSL  ■  mean  sea  level. 
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KM 


MSL,**  ft 


25.0 

29.0 

33.0 

36.9 

40.7 

45.7 


-2.9 

-0.2 

+2.7 

+5.5 

+8.0 

+12.0 


Tide  Gate  Characterization 


19.  The  tide  gate  is  located  on  Back  River.  The  tide  gate  is  con¬ 
trolled  to  open  during  the  flood  tide  and  close  at  high-water  slack.  The 
switching  of  the  tide  gate's  position  is  automatically  controlled  by  sensing 
changes  in  the  water  surface  elevation. 

20.  The  tide  gate  operation  was  simulated  by  defining  Back  River  as  two 
branches:  Branch  2  upestuary  and  Branch  3  downestuary  of  the  tide  gate  (Fig¬ 
ures  2  and  3) .  Flow  through  the  tide  gate  during  the  flood  tide  was  simulated 
as  a  withdrawal  from  Branch  3  and  an  Inflow  into  Branch  2.  During  ebb  tide, 
withdrawals  from  Branch  3  and  inflows  into  Branch  2  were  specified  equal  to 
zero. 

21.  Flow  through  the  tide  gate  was  calculated  assuming  that  the  open 
tide  gate  could  be  modeled  as  a  submerged  weir.  The  bottom  of  the  tide  gate 
structure  is  located  4.145  m  below  MSL.  The  width  of  the  structure  is 
183.5  m.  The  formulae  selected  were  from  King  and  Brater  (1963).  In  refer¬ 
ence  to  Figure  4,  which  is  a  schematization  of  the  tide  gate,  the  flow  Q 
through  the  structure  was  calculated  using  the  relationship: 


where 


Q  -  discharge  through  the  gate,  m  /sec 


Q*  -  discharge  at  the  head  H.  computed  from  the  equation  for  free  dis- 

3  1 

charge,  m  /sec 


Hj  -  head  on  the  upstream  side  of  the  weir,  m 
H.,  -  head  on  the  downstream  side  of  the  weir,  m 
n  ■  1.5  -  exponent  in  the  free  discharge  equation 
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The  free  discharge  equation  used  was 


f'*  -  183.5  C 


(2) 


where 

C  »  2.0  ■  discharge  coefficient 
n  -  1.5 

The  head  heights  and  H^  were  calculated  based  on  the  segment  elevations 

of  Branch  2  and  Branch  3  adjacent  to  the  tide  gate  structure.  The  withdrawals 
from  Branch  3  were  proportioned  to  layers  2  through  5  based  on  their  cross- 
sectional  area. 


Constituent  Reactions 


22.  The  primary  constituent  of  Interest  was  DO.  Salinity  was  simulated 
because  gradients  of  this  constituent  affected  the  hydrodynamics  and  provided 
a  means  to  calibrate  the  model  for  constituent  transport. 

23.  Two  model  formulations  were  simulated.  The  first  was  a  DO  and  BOD 
relationship,  and  the  second  included  the  nutrients  nitrogen  and  phosphorus 
and  the  biological  component  algae. 

24.  Constituent  model  formulations  are  detailed  In  the  report 
"CE-QUAL-W2:  A  Numerical  Two-Dimensional,  Laterally  Averaged  Model  of  Hydro¬ 
dynamics  and  Water  Quality;  User's  Manual"  (Environmental  and  Hydraulics  Labo¬ 
ratories  1986).  The  following  discussion  describes  the  constituent 
interpretation  permitting  specific  application  to  the  Savannah  River  estuary. 
Formulation  1 ;  DO-BOD 

25.  BOD  was  simulated  as  two  components  because  Identified  oxygen 
demanding  substances  varied  greatly  in  their  decay  rates.  BOD  due  to 
effluents  from  the  pulp  and  paper  mill,  secondary  treated  municipal  sources, 
and  the  urban  canal  were  simulated  as  one  component  (BOD) ,  whereas  the  ferrous 
sulfate  (Fe^SO^)  effluent  from  American  Cyanamid  was  simulated  as  a  second 
component  (BOD*).  Both  components  were  modeled  assuming  first-order  decay  but 
with  different  reaction  rates. 

26.  Background  BOD  and  BOD  due  to  pulp  and  paper,  secondary  treated 
municipal,  and  urban  canal  effluents  were  simulated  as  constituent  5,  and 


BOD*  was  simulated  as  constituent  6  In  CE-QUAL-W2.  Constituent  5  is  typically 
called  labile  dissolved  organic  matter  (DOM)  in  CE-QUAL-W2,  and  constituent  6 
is  called  refractory  DOM.  However,  since  the  kinetic  description  of  these  two 
variables  was  Identical  to  those  required  to  simulate  BOD  and  BOD*  in  this 
application,  the  only  modification  required  was  in  the  value  of  the  kinetic 
coefficients  and  in  the  naming  convention.  This  convention  made  use  of  the 
existing  generality  of  CE-QUAL-W2.  The  specification  of  reaction  rate  coeffi¬ 
cients  permitted  the  characterization  of  appropriate  pathways  and  reaction 
rates.  The  pertinent  reaction  rate  coefficients  for  Formulation  1  are  listed 
in  Table  2. 

27.  Other  constituents  simulated  were  temperature,  salinity  (constit¬ 
uent  4),  and  DO  (constituent  12).  The  interfacial  DO  exchange  rate  was  com¬ 
puted  according  to  O'Connor  (1983).  The  transfer  coefficients  were  based  on 
intermediate  (field  average)  scales.  Reaction  rates  were  modified  by  water 
temperature  using  the  formulation  of  Thornton  and  Lessem  (1978). 

Formulation  2 

28.  Formulation  2  included  salinity  (constituent  4),  BOD  (constit¬ 
uent  5) ,  algae  (constituent  7) ,  detritus  (constituent  8) ,  phosphorus  (constit¬ 
uent  9),  ammonia  (constituent  10),  nitrate  +  nitrate  (constituent  11),  and  DO 
(constituent  12).  BOD*  effluent  from  American  Cyanamld  was  simulated  as  a 
separate  component  (constituent  6)  similar  to  Formulation  1. 


Boundary  Conditions 

29.  In  order  to  model  water,  heat,  and  constituent  budgets  accurately 
through  the  simulation  period,  it  was  necessary  to  consider  heat  and  light 
exchange  through  the  water  surface,  water  and  constituent  mass  and  heat 
advected  in  and  out  of  the  estuary,  and  DO  exchange  through  the  water  surface 
and  with  the  bottom  sediments. 

Meteorological  data 

30.  Weather  information  was  acquired  from  NOAA's  Tape  Deck 

No.  1440  WBAN  Hourly  Surface  Observations  for  the  Savannah  Municipal  Airport 
(Station  No.  03822).  Meteorological  data  were  needed  to  compute  equilibrium 
temperature  (ETM) ,  coefficient  of  surface  heat  exchange  (CSHE) ,  solar 
radiation  (SRO),  and  wind  shear  at  the  water  surface.  The  procedure  used  to 
prepare  the  meteorological  data  included:  (a)  extraction  of  only  those  data 


Table  2 

Reaction  Rate  Coefficients  for  Formulation  1 


Parameter 

Card 

Value 

Description* 

TDOMDK 

RI 

0.3 

Maximum  labile  DOM  decay  rate,  1/day  (BOD  decay 
rate) 

TRFRDK 

RJ 

2.0 

Maximum  refractory  DOM  decay  rate,  1/day  (BOD* 
decay  rate) 

TDOMRF 

RJ 

0.0 

Maximum  rate  at  which  labile  DOM  decomposes  to 
refractory  DOM,  1/day  (a  value  of  0.0  was 
specified  since  there  was  no  pathway  between  BOD 
and  BOD*) 

DOMCNT 

RJ 

1.0 

Fraction  of  influent  DOM  that  is  labile** 

D0MT1 

RK 

2.0 

Lower  temperature  bound  of  DOM  decay,  °C 
(BOD  and  BOD*) 

DOMT2 

RK 

20.0 

Lower  temperature  bound  of  maximum  DOM  decay, 

°C  (BOD  and  BOD*) 

D0MK1 

RK 

0.12 

Rate  multiplier  corresponding  to  DOMT1  (BOD  and  BOD*) 

DOMK2 

RR 

0.98 

Rate  multiplier  corresponding  to  DOMT2  (BOD  and  BOD*) 

02D0M 

RO 

1.0 

Oxygen  stoichiometric  equivalent  for  DOM  decay 

(BOD  and  BOD*) 

*  The  first  description  represents  the  general  naming  convention  in 

CE-QUAL-W2,  whereas  the  description  in  parentheses  represents  the  conven¬ 
tion  used  in  this  implementation. 

**  The  code  permits  input  of  DOM  with  subsequent  partitioning  between  labile 
and  refractory  components.  The  application  explicitly  specified  the 
inputs  of  labile  (BOD)  and  refractory  (BOD*)  dissolved  organic  matter; 
thus,  no  additional  partitioning  was  required. 


needed  as  input  to  the  model,  (b)  flagging  missing  data,  (c)  translation 
overpunch  data,  and  (d)  conversion  of  the  hourly  observation  data  to  daily 
averages.  The  computer  codes  described  in  Environmental  Laboratory  (1982) 
were  modified  and  used  for  meteorological  data  preparation. 

31.  The  computations  of  equilibrium  temperature,  coefficient  of  surface 
heat  exchange,  and  solar  radiation  were  based  on  the  Heat  Exchange  Program 
(Program  No.  722-F5-E1010)  available  from  the  US  Army  Engineer  District, 


Baltimore.  Code  modifications  were  necessary  to  ensure  units  consistent  with 
the  model  requirements. 

External  flow  boundary  conditions 

32.  External  flow  boundary  conditions  consisted  of  freshwater  flow  of 
the  Savannah  River  at  the  upestuary  boundary,  two  withdrawals  of  river  water, 
and  8 lx  point  source  discharges.  Each  Inflow  required  a  magnitude,  tempera¬ 
ture,  and  concentration  for  each  modeled  constituent.  The  withdrawals 
required  specification  of  magnitude. 

33.  Freshwater  flow.  The  US  Geological  Survey  (USGS)  publishes  water- 

discharge  records  at  Clyo,  Ga. ,  at  RM  60.9  (USGS  1980),  which  is  26  km 

upstream  from  the  upestuary  boundary.  The  discharge  during  August  1980  was 

3 

rather  uniform  averaging  204  m  /sec  with  a  maximum  and  minimum  of  248  and 
3 

173  m  /sec,  respectively.  The  flow  increase  between  the  gage  and  the  model's 
upestuary  boundary  was  based  on  the  ratio  of  drainage  areas  of  1.072  (LMS 
1983).  Freshwater  flows  were  updated  dally. 

34.  The  dlel  study  provided  freshwater  Inflow  estimates  of  temperature, 
DO,  and  conductivity;  the  center-line  study  results  of  the  upestuary  station 
at  Ebenezer  Landing  provided  freshwater  inflow  estimates  of  the  nutrients  NH^, 
NOj  +  NO^ ,  total  P,  and  BOD^.  Linear  Interpolation  between  sampling  dates 
provided  daily  updates. 

35.  Withdrawals.  Two  withdrawals  of  water  were  simulated  In  the  model 
application:  the  Savannah  Water  Works  (SWW)  Surface  Filtration  Plant  (RM  29, 

Segment  15)  and  the  Beaufort-Jasper  Water  Authority  (BJWA)  Freshwater  Canal 

3 

(RM  39,  Segment  7).  The  SWW  withdrawal  of  1.54  m  /sec  and  the  BJWA  withdrawal 
3 

of  0.26  m  / sec  were  assumed  constant. 

36.  Waste loads.  Point  source  discharges  were  simulated  as  tributaries. 
Three  industries  (American  Cyanamld,  Union  Camp,  and  Continental  Group),  a 
municipal  sewage  treatment  plant  (President  Street  STP),  and  two  urban  canals 
(Springfield  and  Plpemakers)  were  included  as  point  source  loads.  Effluent 
characteristics  were  abstracted  from  LMS  (1983)  and  tabulated  in  Table  3. 

37.  An  ultimate  BOD  to  BOD^  factor  of  1.5  was  assumed  for  the  Preside  t 
Street  STP  and  the  two  urban  canals.  An  ultimate  BOD  to  BOD^  factor  of  2.5 
was  selected  for  the  effluent  from  the  Union  Camp  and  Continental  Group  pulp 
and  paper  mills.  These  estimates  were  abstracted  from  LMS  (1983). 
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Ultimate 

Discharge  BOD  DO  Salinity  Temperature 


38.  American  Cyanamld  did  not  discharge  organic  wastewater,  but  the 

f  j 

chemical  manufacturer  did  discharge  ferrous  sulfate.  The  oxidation  of  Fe  to 

j  (  | 

Fe  occurs  rapidly  and  depletes  the  DO  of  the  receiving  waters.  LMS  (1983) 
simulated  the  American  Cyanamld  iron  wasteload  as  an  immediate  oxygen  demand 
using  a  DO  deficit  in  their  model.  In  contrast,  the  present  study  simulated 
ferrous  sulfate  as  a  separate  component,  undergoing  first-order  decay  requir¬ 
ing  DO  similar  to  BOD,  but  at  an  increased  rate. 

External  head  boundary  elevation 

39.  A  Fisher  and  Porter  Type  150  surface  elevation  level  recorder  was 
located  at  the  downstream  study  boundary  at  RM  5.9.  The  gage  data  were  used 
to  specify  the  downstream  external  head  boundary  elevation.  Analysis  of  the 
tide  gage  data  consisted  of  data  editing  to  fill  in  short  sequences  of  missing 
data,  filtering  to  remove  high-  and  low-frequency  trends  from  the  data,  and 
harmonic  analysis  for  tidal  constituents. 

40.  Data  editing.  For  single  missing  data  values,  an  average  of  adja¬ 
cent  tidal  heights  was  substituted.  For  longer  intervals  of  missing  data,  the 
preceding  tidal  cycle  water  surface  elevations  were  superimposed  over  the 
missing  interval,  and  tidal  heights  were  visually  estimated.  The  manual  prep¬ 
aration  of  tide  gage  data  resulted  in  a  tidal  record  extending  from  24  July 
1980  at  1400  EST  through  27  August  1980  at  0654  EST  with  8,100  observations. 

41.  Filtering.  A  digital  band-pass  filter  was  applied  to  attenuate 
high  and  low  frequencies.  The  period  range  considered  in  the  harmonic  analy¬ 
sis  was  approximately  12  to  24  hr.  Because  of  the  large  tidal  range,  absence 
of  major  passing  synoptic  weather  fronts,  and  uniform  freshwater  flow  during 
the  recording  period,  filtering  was  not  mandatory.  However,  filtering  did 
result  in  smaller  deviations  between  measured  and  predicted  elevations  because 
noise  not  accounted  for  in  the  descriptive  equation  for  water  surface  eleva¬ 
tion  was  attenuated.  An  eight-pole  Butterworth  filter,  characterized  by  a 
smooth  power  gain  with  maximum  flatness  in  the  passband  and  the  stop  band 
along  with  a  reasonably  sharp  cutoff,  was  used  (Outlaw  and  Butler  1982). 
Half-power  frequencies  were  4.63  *  10  ^  Hz  and  7.72  *  10*"^  Hz  or  6  hr  and 

36  hr,  respectively. 

42.  Harmonic  analysis.  The  arithmetic  mean  of  the  tide  record  was  sub¬ 


tracted  from  each  datum,  and  the  result  was  interpreted  as  a  deviation  from 
MSL,  the  datum  used  in  the  Savannah  River  Estuary  study. 


43.  The  elevation  of  the  Savannah  River  Estuary  tide  was  represented  by 


h(t)  “  ^  Aj  co^  cj+  <f>  )  (3) 

j-1 


where 


h 

t 

m 


A 

u 

♦ 


3 

j 

j 


elevation  at  time  t 

time  from  some  initial  time  (1400,  24  July  1980,  in  the  study) 

number  of  tidal  constituents 

amplitude  of  the  j**1  tidal  constituent 

angular  frequency  of  the  j  tidal  constituent 

phase  of  the  j**1  tidal  constituent 


An  observed  surface  elevation  datum  h^(t)  can  be  represented  by 


ht(t) 


2 

j-1 


(a^  cos 


V 


+  bj  sin 


“jt) 


ei(t) 


(4) 


where 

Sj ,  bj  -  coefficients 
e^(t)  -  error  in  observed  datum  i 


For  each  tidal  constituent,  the  amplitude  A^  and  the  phase 
determined  from 


can  be 


A 


j 


i 


(5) 

(6) 


44.  A  least  squares  analysis  was  conducted  to  estimate  the  coefficients 
a^  and  ,  which  minimized  the  sum  of  squares  of  the  N  -  8,100  error 
terms 


(v-i) i  v 

x  >  i-i 


The  N^f  Sj,  K^,  0^,  P , ,  and  tidal  harmonics  were  estimated.  The  tidal 

harmonic  constituents  represent  selected  effects  of  the  gravitational  effects 
of  the  earth,  moon,  and  sun.  The  subscript  1  represents  diurnal  periods,  and 
the  subscript  2  represents  semidiurnal  periods. 

Result 

45.  Component  amplitudes  and  phases  for  the  tidal  record  are  tabulated 
in  Table  4.  The  observed,  filtered,  and  simulated  tidal  elevations  are  dis¬ 
played  in  Figure  5.  Equation  3,  using  the  tidal  constituent  values  in 
Table  4,  was  used  in  the  subroutine  TVDS  to  provide  the  downstream  external 
head  boundary  elevation  (ELR)  for  Branch  1. 

External  head  boundary  profiles 

46.  On  6,  13,  19,  and  26  August  1980,  surface-to-bottom  profiles  of 
temperature,  DO,  and  conductivity  were  made  at  1-m  intervals  at  the  downstream 
boundary.  Uater  samples  for  NH^,  NO2  +  NO^,  total  P,  and  BOD^  analyses  were 

Table  4 

Calculated  Tidal  Components 


Tidal 

Component 


Reference 

Time 

Period,  h 

Amplitude,  m 

Phase,  rad 

24  July  1980 

12.421 

1.030 

1.892 

1400  EST 

12.658 

0.205 

-0.061 

12.000 

0.228 

-2.533 

23.934 

0.046 

0.497 

25.819 

0.080 

1.073 

24.066 

0.115 

1.437 

11.967 

0.099 

1.100 
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Figure  5.  Prototype,  filtered,  and  simulated  water 
surface  elevation 

collected  at  four  depths:  l  m  below  the  surface,  1  m  above  the  bottom,  and 
two  depths* intermediate  to  surface  and  bottom.  Profiles  were  made,  and  sam¬ 
ples  were  collected  every  hour  for  13  consecutive  hours. 

47.  Interest  was  in  both  intratidal  and  intertidal  variations  over  the 
month.  The  desired  external  head  boundary  profiles  required  intratidal 
updates  extending  over  the  entire  month  and  at  a  vertical  resolution  corre¬ 
sponding  to  the  model  layers.  Models  specifying  external  head  boundary 
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concentrations  were  needed  to  specify  concentrations  for  times  and  dates  for 
which  no  data  existed. 

48.  Preparation  of  the  profile  data  consisted  of  conversion  of  conduc¬ 
tivity  to  salinity,  relating  sample  depth  to  MSL,  interpolation  over  model 
grid  layers,  and  harmonic  analysis. 

49.  Conversion  of  conductivity  to  salinity.  Conversion  of  conductivity 
(mmho/cm)  to  salinity  (ppt)  was  achieved  by  using  a  piecewise  linear  relation¬ 
ship  abstracted  from  a  nomogram  published  by  Hydrolab  (undated) .  The  rela¬ 
tionship  used  was: 

Conductivity,  mmho/cm  _ Salinity _ 


COND 

< 

10 

SAL 

m 

0. 55*C0ND 

10 

£ 

COND 

< 

20 

SAL 

m 

0. 62*COND-0. 7 

20 

< 

COND 

< 

30 

SAL 

m 

0.67*COND-1.7 

30 

£ 

COND 

< 

40 

SAL 

m 

0.69*COND-2.3 

40 

COND 

< 

50 

SAL 

m 

0. 75*COND-4. 7 

50 

< 

COND 

SAL 

m 

0. 74*COND-4. 2 

50.  Relation  of  sample  depth  to  MSL.  Samples  were  collected  at  depths 
relative  to  the  water  surface.  However,  the  water  surface  elevation  changed 
within  the  13-hr  sampling  period  as  well  as  between  periods.  Equation  3  with 
tidal  constituents  for  the  tide  gauge  at  RM  5.9  was  used  to  specify  the  water 
surface  elevation  at  the  time  of  sample  collection  and  permitted  specification 
of  sample  depth  relative  to  MSL. 

51.  Interpolation  over  model  grid.  Linear  interpolation  was  used  to 
assign  a  constituent  concentration  to  each  of  the  model  grid  layers.  Grid 
layers  deeper  than  measured  data  were  assigned  values  equal  to  the  adjacent 
Interior  measured  data. 

52.  Harmonic  analysis.  Examination  of  the  data  revealed  periodicities 
for  all  constituents  and  a  linear  decrease  with  time  for  temperature  between 
dates.  The  average  tidal  period  for  August  1980  abstracted  from  (NOAA)  Tide 
Tables  (NOAA  1980)  for  Savannah  was  12.42  hr,  the  period  of  the  M2  tidal 
constituent.  Therefore,  least  squares  were  used  to  fit  a  sinusoidal  with  a 
period  of  12.42  hr  (angular  frequency  “  0.5059  rad/hr)  to  the  observed  data 
for  each  layer.  A  linear  regression  term  was  Included  for  temperature. 

53.  The  model  used  for  temperature  was 

T  -  T  +  Bt  +  A  cos  (wt  -  0)  (8) 

o 
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For  salinity,  DO,  BOB,.,  NH^,  NO,, ,  and  F,  the  model  used  was 

C  ■  C  +  A  cos  (w t  -  0)  (9) 

o 

where 

T  -  temperature,  °C 
Tq  ■  Intercept,  °C 
8  *  regression  coefficient,  °C/hr 
t  ■  time,  hr 

C  -  constituent  concentration,  ppt  or  g/m9 
Cq  ■  constituent  mean  concentration,  ppt  or  g/m9 
A  -  amplitude,  °C  or  g/m9 
a)  ■  angular  frequency,  rad/hr 
0  *  phase,  rad 

54.  The  Naval  Surface  Weapons  Center  (NSWC)  (1981)  least  squares  sub¬ 
routines  were  used  to  fit  the  selected  models. 

55.  In  Table  5  are  listed  the  Intercept,  regression  coefficient,  ampli¬ 
tude,  phase,  and  R2  value  of  the  temperature  model  for  each  layer.  The 

R2  ,  or  coefficient  of  determination,  measures  the  portion  of  the  variation  in 
the  temperature  that  is  attributed  to  the  model  rather  than  to  random  error. 

R2  is  defined  as  follows: 


(10) 


where 

T^  *  observed  temperature  datum 

A 

T^  *  predicted  temperature  based  on  model 
T  ■  mean  of  temperature  data 

56.  In  Tables  6  through  11  are  listed  the  mean,  amplitude,  phase,  and 
R2  value  for  salinity,  DO,  BODj,  NH^-N,  N02  +  N03~N,  and  total-P. 

57.  Examination  of  Tables  5  through  11  reveals  that  temperature  and 
salinity  were  adequately  characterized  by  the  sinusoidal  function,  that  the 
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Table  5 

Intercept,  Regression  Coefficient,  Amplitude,  Phase,  and  Coefficient  of 


2 

De termlna t ion  (R  )  for  the  Temperature  Profile  Head  Boundary  Condition 


Intercept 

Regression  Coefficient 
_2 

Amplitude 

Phase* 

1 

2  i 

1  Layer 

°C 

10  °C/hr 

°C 

rad 

R 

m 


Table  6 

2 

Mean,  Amplitude,  Phase,  and  Coefficient  of  Determination  (R  )  for 


the  Salinity  Profile  Head  Boundary  Condition 


Amplitude 
t 


Phase* 

rad 


Table  7 

2 

Mean,  Amplitude,  Phase,  and  Coefficient  of  Determination  (R  )  for 


the  DO  Profile  Head  Boundary  Condition 


Amplitude 


Layer 


Phase* 

rad 


*  Reference  time  -  1  August  1980,  0000  EST. 


Table  8 

Kean,  Amplitude,  Phase,  and  Coefficient  of  Determination  (RZ)  for 
the  BOD j  Profile  Head  Boundary  Condition 


Mean 


Amplitude 


Phase* 

rad  RZ 


Layer 


Table  9 


ll 


Mean,  Amplitude,  Phase,  and  Coefficient  of  Determination  (R  )  for 


Table  10 


2 

Mean,  Amplitude,  Phase,  and  Coefficient  of  Determination  (R  )  for 
the  NO2  +  NO^-N  Profile  Head  Boundary  Condition 


Mean 

g/m3 


Amplitude 

g/m3 


Phase* 

rad 


1 

! 

( 


I 


Layer 


Table  11 

2 

Mean,  Amplitude,  Phase,  and  Coefficient  of  Determination  (R  )  for 
the  Total-P  Profile  Head  Boundary  Condition 


Mean 

-J 

Amplitude 

9 

Phase* 

Layer 

g/m 

g/m-1 

rad 

R2 

2 

0.07 

0.01 

-3.00 

0.27 

3 

0.07 

0.01 

-2.76 

0.11 

4 

0.07 

0.01 

-2.60 

0.05 

5 

0.08 

0.00 

-1.86 

0.02 

6 

0.08 

0.01 

-1.26 

0.05 

7 

0.08 

0.01 

-1.13 

0.10 

8 

0.09 

0.02 

-1.04 

0.13 

9 

0.10 

0.02 

-0.89 

0.16 

10 

0.11 

0.03 

-0.94 

0.21 

11 

0.12 

0.05 

-1.04 

0.21 

12 

0.14 

0.06 

-1.11 

0.19 

13 

0.16 

0.07 

-1.16 

0.18 

14 

0.17 

0.07 

-1.28 

0.17 

*  Reference  time  ■  1  August  1980,  0000  EST 


sinusoidal  function  accounts  for  45  percent  of  the  variability  of  DO  and 
NO2  +  NO^,  and  that  the  selected  sinusoidal  function  poorly  fits  the  BOD,  NH^, 
and  P  data.  In  Figures  6  through  25,  prototype  and  simulated  data  values  for 
temperature,  salinity,  BOD^,  and  DO  are  displayed. 

Sediment  oxygen  demand  (SOD) 

58.  SOD  demand  rates  were  determined  at  RM  5.9  (Segment  33),  RM  16.3 

(Segment  25),  and  RM  21.5  (Segment  21)  from  1  through  5  October  1980 

(Table  12) .  The  results  were  linearly  Interpolated  over  the  model  segments 
2 

assuming  1.2  g/m  / day  for  Segment  2. 


Table  12 


Sediment 


gen  Demand  (SOD) 


Segment 

21 


Mean  Rate,  0, 

2  i 
g/m  /day 


Initial  Conditions 

59.  Initial  conditions  consist  of  water  surface  elevation,  temperature, 
and  constituent  concentrations  corresponding  to  the  simulation  starting  time. 
The  simultaneous  solution  of  the  combined  momentum  and  continuity  equations 
establishes  a  consistent  flow  field  within  a  few  wave  travel  times  down  the 
estuary;  however,  the  establishment  of  a  conservative  constituent  field  con¬ 
sistent  with  the  applied  boundary  conditions  requires  an  amount  of  time 
approaching  the  residence  time  (Buchack  and  Edinger  1982). 

60.  Initial  conditions  must  be  specified  for  each  ^odel  cell.  The  pro¬ 
cedure  used  was  to  Initialize  all  cells  to  one  constituent  value,  simulate  for 
a  period  of  time  sufficient  for  the  constituent  field  to  equilibrate  to  the 
boundary  conditions,  and  use  the  resulting  simulated  cell  constituent  values 
as  initial  conditions  for  subsequent  runs.  An  extended  simulation  with  moni¬ 
toring  of  transient  cell  constituent  concentrations  provided  an  estimate  of 
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Figure  6.  Measured  temperatures*  August  1980 
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Figure  7.  Measured  and  calculated  temperatures* 
6  August  1980 
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Figure  8. 
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Figure  9.  Measured  and  calculated  temperatures 
19  August  1980 
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Figure  12.  Measured  and  calculated  salinity 
6  August  1980 
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Figure  13.  Measured  and  calculated  salinity 
13  August  1980 
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Figure  18.  Measured  and  calculated  BOD,, 
13  August  1980 
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Figure  19.  Measured  and  calculated  BOD,, 
19  August  1980  J 
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Figure  26.  Transient  response,  simulated  tldally 
averaged  salinity,  Segment  21 

In  Segments  21,  27,  and  29  at  the  surface,  middepth,  and  bottom.  Equilibrium 
of  salinity  was  achieved  In  10  to  20  tidal  cycles.  Figures  29  through  31 
represent  the  tldally  averaged  salinity  profile  after  1,  25,  and  50  tidal 
cycles.  Figures  32  through  34  present  the  corresponding  tldally  averaged 
velocity.  For  the  Savannah  River  Estuary  simulations,  a  restart  file  based  on 
28  days  (54  tidal  cycles)  was  used.  Test  runs  were  initiated  from  the  restart 
file,  and  these  ran  for  5  days  (10  tidal  cycles)  before  comparisons  were  made. 
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Figure  28.  Transient  response,  simulated  tidally 
averaged  salinity.  Segment  29 
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9.  Simulated  tidally  averaged.  Figure  30.  Simulated  t 

1  salinity  profile.  Segment  21  vertical  salinity  prof 


PART  III:  RESULTS 


Hydrodynamic  Adjustment 

i 

t 

( 

Tides  and  velocities  1 

“■  ^  “  11  I 

Cl.  Two  results  of  the  simulations  were  used  to  achieve  hydrodynamic  j 

adjustment:  spatial  and  temporal  variation  of  tide  height  and  the  pattern  of 

t 

circulation  throughout  the  tidal  cycle.  The  numerical  model  provided  one  ‘ 

coefficient,  the  Chezy  coefficient,  that  had  detectable  influence  on  tidal 
variation  and  velocity  fields  and  could  be  varied  to  achieve  hydrodynamic 
adjustment.  The  Chezy  coefficient  relates  vertical  velocity  gradients  to  bot¬ 
tom  and  interlayer  stresses. 

63.  Initial  simulations  revealed  that  simulated  tidal  amplitudes  were 
nearly  equal  to  measured  (surface  elevation  recorder  at  RM  11,  Segment  29)  and 
predicted  (NOAA  Tide  Tables  (1980))  amplitudes,  but  that  tidal  phase  lags 
occurred  and  increased  upestuary  (Table  13).  The  phase  errors  varied  with  the 
Chezy  coefficient,  but  phase  errors  could  not  be  eliminated  by  manipulation  of 
the  Chezy  coefficient  alone. 

64.  Simulated  velocities  increased  with  increased  values  for  the  Chezy 
coefficient,  but  simulated  surface  velocities  were  consistently  less  than 
those  measured.  Simulation  results  also  indicated  that  tidal  fluctuations 
increased  in  magnitude  at  the  upstream  boundary  with  the  increased  Chezy 
coefficient.  The  upstream  boundary  at  Ebenezer  Landing  is  upestuary  of  the 
fall  line,  and  thus  the  water  surface  elevation  should  not  tidally  fluctuate 
in  the  simulations.  In  order  to  eliminate  simulated  tidal  fluctuation  at  the 
upstream  boundary  and  to  minimize  tidal  phase  errors,  the  coding  was  modified 
to  make  the  Chezy  coefficient  equal  to  30  in  the  navigation  channel  and  sedi¬ 
ment  basin  and  25  in  the  nondredged  channel. 

65.  The  observation  that  simulated  tidal  elevations  were  nearly  equiva¬ 
lent  to  measured  values  whereas  simulated  tidal  velocities  were  less  than  mea¬ 
sured  values  implied  that  the  system  geometry  was  in  error.  The  detailed 
"Savannah  Harbor,  Georgia,  Annual  Survey"  (SAS  1981)  and  the  detailed  "Naviga¬ 
tion  Charts,  Savannah  River,  Georgia  and  South  Carolina"  (SAS  1980)  that  were 
used  to  specify  geometry  for  Branches  1  and  3  were  considered  adequate.  How¬ 
ever,  the  Nautical  Chart  11514  (NOAA  1978)  that  was  used  to  specify  Branch  2 
geometry  contained  little  detail  for  the  non-navigable  channel.  The  influence 


of  underestimating  the  storage  of  Branch  2  was  examined  by  systematically 
increasing  the  section  widths  and  noting  the  resulting  ebb  velocities  in  Front 


River.  Increasing  Branch  2  storage  by  20  percent  increased  ebb  velocities  by 
less  than  4  percent.  Because  of  the  apparent  Insensitivity  of  resulting  ebb 
velocities  to  moderate  increases  in  storage,  the  estuarine  geometry  that  could 
be  abstracted  from  the  harbor  survey,  navigation,  and  nautical  charts  (SAS 
1980  and  1981;  NOAA  1977  and  1978)  was  used  without  modification.  Subsequent 
investigations  revealed  that  extensive  storage  occurs  through  overbank 
flooding.  Overbank  flooding  was  not  included  in  the  model. 

66.  Simulated  and  measured  tidal  elevations  at  RM  11  (Segment  29)  are 
displayed  in  Figure  35.  Simulated  and  measured  surface,  middepths,  and  bottom 
velocities  at  RM  11.4  (Segment  29)  and  RM  17.2  (Segment  24)  are  illustrated  in 
Figures  36  through  41.  Examination  of  Figure  35  reveals  that  tidal  amplitude 
and  phase  were  reproduced  in  downestuary  segments.  However,  as  discussed 
previously,  tidal  phase  errors  occurred  at  upestuary  locations  (Table  13). 
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Figure  35.  Simulated  and  measured  tidal  elevations  at 

Segment  29 
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Figure  36.  Simulated  and  measured  surface  velocities  at  Segment  29 
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Figure  37.  Simulated  and  measured  middepth  velocities  at  Segment  29 
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Figure  38.  Simulated  and  measured  bottom  velocities  at  Segment  29 
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Figure  40.  Simulated  and  measured  middepth  velocities  at  Segment  24 
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Comparison  of  simulated  and  measured  velocities  indicates  greater  differences 
at  Segment  29  than  at  Segment  24,  greater  differences  during  ebb  (velocities 
positive)  than  during  flood,  and  greater  differences  at  the  surface  than  at 
the  bottom. 

67.  The  distinction  between  right  and  left  bank  simulated  velocities  in 
Figures  36  through  41  was  due  to  the  fact  that  simulated  and  measured  veloci¬ 
ties  were  matched  with  respect  to  time  and  depth  of  measurement.  The  time 
required  to  sample  both  the  left  and  right  bank  was  approximately  30  min;  the 
comparison  of  measured  and  simulated  velocities  corresponds  to  within  3  min, 
with  the  result  that  left  and  right  bank  comparisons  may  differ  by  30  min. 

The  sampling  was  conducted  adjacent  to  the  navigation  channel,  and  the  mea¬ 
sured  depths  differed  on  either  side  of  the  channel;  therefore,  bottom  and 
middepth  measurements  differed  in  elevation  and  frequently  corresponded  to 
different  model  layers. 

68.  Variations  between  observed  and  predicted  tidal  velocities  and 
phases  were  attributed  in  part  to  Inaccuracies  in  the  model  geometry.  Two 
observations  need  to  be  considered  in  future  estuarine  applications.  First, 
the  navigation  charts  used  for  specification  of  the  system  geometry  typically 
recorded  depths  no  shallower  than  MLW.  The  specification  of  segment  widths  at 
elevations  no  higher  than  MLW  is  duplicated  by  the  computer  code  GEDA  at 
higher  elevations.  This  results  in  an  underestimate  of  cross-sectional  areas 
during  the  tidal  cycle  at  elevations  greater  than  MLW.  Topographic  maps  did 
not  provide  the  resolution  required  for  the  Savannah  River  application. 

Second,  extensive  areas  of  the  Savannah  River  Estuary  are  subject  to  overbank 


flooding  during  the  tidal  cycle,  and  nautical  charts  may  be  in  error  for  non- 
navigable  streams.*  These  observations  suggest  that  recent  detailed  topo¬ 
graphic  surveys  be  used  or  that  personal  observation  be  employed  to  define  the 
areal  extent  and  depth  of  overbank  flooding  in  order  to  accurately  reproduce 
the  systems  geometry. 

Salinity 

69.  Figure  42  illustrates  a  comparison  of  simulated  and  measured  salin¬ 
ity  profiles  for  Segments  24,  27,  and  29  for  four  sampling  dates.  The  mea¬ 
sured  salinities  were  abstracted  from  the  center-line  water  quality  survey 


*  Personal  Communication,  9  August  1985,  Billy  H.  Johnson,  Research  Hydraulic 
Engineer,  US  Army  Engineer  Waterways  Experiment  Station,  Vicksburg,  Miss. 
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Figure  42.  Simulated  and  measured  salinities  at 
Segments  24,  27,  and  29 
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(Pennington  and  Bond  1981).  Figures  43  through  46  represent  snapshots  of  the 
simulated  velocity  structure  and  water  surface  elevation  corresponding  to  the 
time  of  profile  measurement.  The  simulated  salinities  closely  approximate  the 

measured  values  on  August  5  and  12  but  deviate  on  August  20  and  27.  Adjust- 

2 

ment  of  the  horizontal  dispersion  coefficient  between  1  and  10  m  /sec  had  no 
detectable  effect  on  salinity  concentrations. 

70.  Two  possible  reasons  are  postulated  to  account  for  these  apparent 
deviations.  The  first  is  that  the  center-line  surveys  required  4  hr  for  their 
completion.  The  simulation  snapshots  for  profile  comparison  were  made  at  the 
arithmetic  average  of  sample  collection  times  resulting  in  possible  deviations 
of  up  to  2  hr  between  simulated  and  measured  values.  The  Savannah  River  Estu¬ 
ary  is  quite  dynamic  with  observed  salinity  variations  at  some  stations 
exceeding  S  ppt  over  a  2-hr  period.  The  second  reason  is  the  possible  exis¬ 
tence  of  errors  in  the  advectlve  transport  because  of  unaccounted  overbank 
storage  or  discrepancies  in  geometry.  The  discrepancy  between  measured  and 
simulated  velocities  supports  this  hypothesis. 


71.  The  model  test  program  consisted  of  two  levels  of  constituent  for¬ 
mulations.  The  first  level  consisted  of  only  DO  and  BOD  relationships, 
whereas  the  second  included  the  addition  of  the  nutrients  nitrogen  and 
phosphorus  and  the  biological  component  algae.  Experiments  were  conducted 
using  the  DO-BOD  formulation  to  demonstrate  the  utility  of  CE-QUAL-W2  for 
investigating  changes  in  upstream  regulation  and  evaluating  permit  applica¬ 
tions.  Sensitivity  analyses  were  conducted  using  both  levels  of  kinetic 
formulations  in  order  to  investigate  the  importance  of  tidal  boundary  constit¬ 
uent  concentrations. 

Formulat ion  1 :  DO-BOD 

72.  Calibration.  Examination  of  Figures  47  through  49  reveals  that 
simulated  and  measured  profiles  of  temperature,  BOD,  and  DO  are  nearly  equiva¬ 
lent.  The  results  were  quite  good  considering  the  difficulty  in  specifying 
head  boundary  constituent  profiles. 

73.  Head  boundary  sensitivity  experiment.  The  head  boundary  constit¬ 
uent  profiles  were  based  on  the  1980  survey  results.  A  factorial  sensitivity 
experiment  was  conducted  in  order  to  evaluate  the  effects  of  imprecise 
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Figure  45.  Simulated  velocity  structure  and  water  surface 
elevation,  10:30  AM,  20  August  1980 
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specification  of  DO  and  BOD  head  boundary  concentrations  on  the  simulated 
tidally  averaged  DO  concentration  of  the  whole  modeled  system. 

74.  The  results  of  the  experiment  are  tabulated  in  Table  14.  Examina¬ 
tion  of  Table  14  reveals  no  interaction  between  head  boundary  concentrations 
of  DO  and  BOD  using  the  response  variable  of  tidally  averaged  DO.  Increasing 
the  head  boundary  BOD  concentration  100  percent  decreased  the  tidally  averaged 
DO  concentration  1  percent;  decreasing  the  head  boundary  DO  concentration  50 
percent  decreased  the  tidally  averaged  DO  concentration  10  percent.  Experi¬ 
ments  discussed  in  the  next  section  of  this  report  revealed  that  a  25-percent 
decrease  in  freshwater  flow  and  a  50-percent  increase  in  BOD  loads  had  minor 
effects  on  downstream  DO  and  concentrations.  The  results  of  the  sensitivity 
experiments  coupled  with  the  results  of  the  decreased  flow  and  increased  BOD 
waste  load  experiments  imply  that  significant  errors  were  not  introduced  by 
not  modifying  the  head  boundary  concentrations  of  DO  and  BOD. 

75.  Decreased  flows  and  increased  waste  loads.  Figure  50  displays  con¬ 
tour  plots  of  tidally  averaged  DO  concentrations  under  four  simulated  experi- 

3 

mental  conditions:  (a)  minimum  release  of  173  m  /sec  (6,100  cfs)  from  New 
Savannah  Bluff  Dam,  (b)  25-percent  decrease  in  freshwater  flow,  (c)  50-percent 
increase  in  waste  loads,  and  (d)  a  combination  of  a  25-percent  decrease  in 
freshwater  flow  and  a  50-percent  increase  in  waste  loads.  A  50-percent 
increase  in  waste  loads  was  simulated  by  injecting  126  g/sec  BOD  between 
RM  21.1  (the  US  Highway  Houlihan  Bridge)  and  RM  24.5.  Examination  of  Fig¬ 
ure  50  reveals  a  DO  minimum  near  the  downstream  boundary  under  all  four  simu¬ 
lated  experimental  conditions.  Decreased  flows  resulted  in  a  slight  upstream 
migration  of  the  DO  contours,  and  increased  waste  discharges  resulted  in  a 
slight  decrease  of  the  DO  in  the  bottom  waters  near  Segment  20.  A  combination 
of  decreased  flows  and  increased  waste  loads  resulted  in  conditions  similar  to 
decreased  flows  alone.  The  results  of  the  experiment  suggest  that  the  DO 
concentration  of  the  Savannah  River  is  rather  insensitive  to  moderate 
decreases  in  flow  or  increases  in  waste  discharge. 

Formulation  2:  DO,  algae,  and  nutrients 

76.  Formulation  2  was  an  expansion  of  Formulation  1  that  included  algae 
(constituent  7),  detritus  (constituent  8),  phosphorus  (constituent  9),  ammonia 
(constituent  10),  and  nitrite  +  nitrate  (constituent  11).  The  objective  of 
the  formulation  was  to  evaluate  the  relative  model  sensitivity  of  head  bound¬ 
ary  algal  concentrations  in  comparison  with  the  sensitivity  of  the  maximum 
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Table  14 

Relationship  of  Tidally  Averaged  DO  Concentration  to  Head 
Boundary  Profile  Concentrations  of  DO  and  BOD  Demand 


*  DO  and  BOD  head  boundary  profile  concentration  expressed  as  a  multiple  of 
measured  concentrations. 


TRNOARO  FLOW  EXISTING  WASTE  LOROS  0. ?5»STRNORRO  FLOW  EXISTING  WASTE  LOROS 


algal  growth  rate.  Each  measured  nutrient  was  assumed  to  be  in  a  form  avail¬ 


able  to  the  algae  and  consistent  with  model  formulations.  Reaction  rate 
coefficients  were  selected  equal  to  the  midrange  values  tabulated  in  the 
report  "CE-QUAL-W2:  A  Numerical  Two-Dimensional,  Laterally  Averaged  Model  of 
Hydrodynamics  and  Water  Quality;  User's  Manual"  (Environmental  and  Hydraulics 
Laboratories  1986) . 

77.  Algal  and  detrltal  boundary  conditions  were  not  known.  The  algal 

and  detrltal  flow  boundary  conditions  were  each  set  equal  to  0.5  g  dry 

3  3 

weight /m  ;  the  detrltal  head  boundary  concentration  was  set  equal  to  0.5  g/m  . 

For  comparison,  Thomann  and  Fitzpatrick  (1982)  used  a  phytoplankton  concentra- 

3  3 

tion  of  0.7  g  carbon/  m  (1.6  g  dry  weight/m  )  at  the  upstream  boundary  (Cain 

3  3 

Bridge)  and  0.5  g  carbon/m  (1.1  g  dry  welght/m  )  at  the  downstream  boundary 
(Chesapeake  Bay)  of  their  Potomac  estuary  application. 

78.  The  sensitivity  analysis  consisted  of  a  factorial  experiment  con¬ 
ducted  at  three  levels  of  algal  boundary  concentrations  (1.25,  2.5,  and 

3 

5.0  g/m  )  and  three  levels  of  algal  maximum  growth  rate  (0.25,  0.5,  and 
1.0  e/day) .  The  response  variable  was  defined  as  the  tidally  averaged  algal 
concentration  (see  Table  15) . 

Table  15 

Relationship  of  Tidally  Averaged  Algal  Concentration  to  Head  Boundary 
Algal  Profile  Concentration  and  Algal  Maximum  Growth  Rate 

Algal  Head  Boundary  Maximum  Algal  Growth  Tidally  Averaged  Algal 

3  3 

Concentration,  g/m  _ Rate,  1/day  Concentration,  g/m 


79.  Examination  of  the  results  tabulated  in  Table  15  reveals  no  inter¬ 
action  between  head  boundary  algal  concentration  and  algal  maximum  growth 
rate.  Tidally  averaged  algal  concentration  was  linearly  related  to  both  head 
boundary  concentration  and  maximum  algal  growth  rate.  A  100  percent  error  in 
estimation  of  algal  head  concentration  would  result  in  an  average  algal  mass 
error  of  65  percent,  whereas  a  100  percent  error  in  estimation  of  maximum 
algal  growth  rate  would  result  in  an  average  algal  mass  error  of  11  percent. 
Light  limitation  consistently  was  the  factor  predicted  to  limit  algal  growth. 


PART  IV:  CONCLUSIONS  AND  RECOMMENDATIONS 


Conclusions 


General 

80.  The  2-D  laterally  averaged  water  quality  code  CE-QUAL-W2  character¬ 
ized  the  observed  longitudinal  and  vertical  gradients  of  the  Savannah  River 
Estuary. 

Water  quality  data  resources 

81.  The  data  base  was  adequate  for  model  calibration  of  the  DO-BOD 
formulation  but  inadequate  for  model  confirmation.  A  second,  independent  data 
set  for  model  confirmation  did  not  exist. 

82.  The  absence  of  algal,  detrital  (or  particulate,  nonliving,  organic 
matter) ,  and  inorganic  phosphorus  data  precluded  valid  simulation  of  the 
algae-nutrient-dissolved  oxygen  formulation. 

83.  The  application  provided  the  following  general  comments  concerning 
the  existing  water  quality  data  base: 


a.  The  center-line  survey,  consisting  of  16  center-line  stations 
sampled  on  four  dates,  provided  the  data  used  for  vertical  pro¬ 
file  comparisons.  The  number  of  center-line  stations  on  the 
main  channel  (Savannah  River,  Front  River,  and  North  Channel) 
was  more  than  adequate.  The  sampling  effort  could  have  been 
more  effectively  distributed  with  six  center-line  main  channel 
stations  extending  from  Ebenezer  Landing  downestuary  to  the 
mouth  of  the  Savannah  River  at  Fort  Pulaski  and  two  center-line 
stations  on  the  Back  River.  One  center-line  site  above  RM  25 
would  be  adequate. 

b.  The  four  sediment  oxygen  demand  study  sites  were  limited  in 
spatial  extent.  The  study  sites  need  to  be  extended  to  corre¬ 
spond  with  the  center-line  stations  on  the  Savannah  River, 

Front  River,  Back  River,  and  the  North  Channel. 

Geometric  schematlzation 

84.  Acquisition  of  cross-sectional  information  was  a  critical  problem. 
Navigation  and  nautical  charts  infrequently  recorded  bottom  elevations.  The 
specification  of  segment  widths  at  elevations  no  higher  than  MLW  is  duplicated 
by  the  computer  code  GEDA  at  higher  elevations.  This  results  in  an  underesti¬ 
mate  of  cross-sectional  areas  during  the  tidal  cycle  at  surface  elevations 
greater  than  MLW.  Topographic  maps  did  not  provide  the  resolution  required. 

85.  Extensive  areas  of  the  Savannah  River  Estuary  are  subject  to  over¬ 
bank  flooding  during  the  tidal  cycle.  The  depth  and  extent  of  overbank 


flooding  can  be  only  tangentially  Inferred  from  the  nautical  charts  based  on 
shading  representing  vegetatlonal  cover. 

86.  Instantaneous  water  surface  elevations  extending  over  several 
layers  resulted  in  numerical  instabilities.  Therefore,  a  surface  layer  thick¬ 
ness  of  5  m  was  selected  with  remaining  layers  1  m  thick. 

Tide  gate  characterization 

87.  Tide  gate  operation  was  successfully  simulated  by  defining  Back 
River  as  two  branches:  one  upestuary  and  one  downestuary  of  the  tide  gate. 
Flow  through  the  open  tide  gate  was  modeled  as  a  submerged  weir. 

Constituent  reactions 

88.  The  existing  generality  of  the  water  quality  constituent  formula¬ 
tions  in  CE-QUAL-W2  permits  user  latitude  in  selecting  constitutents  to  simu¬ 
late  and  in  specifying  their  interactions.  The  generality  was  reflected  in 
the  concurrent  simulation  of  BOD  and  BOD*  (Fe^SO^) . 

Boundary  conditions 

89.  The  dlel  study  conducted  at  Ebenezer  Landing  that  was  used  to  spec¬ 
ify  temperature,  DO,  and  salinity  at  the  upstream  flow  boundary  revealed 
little  diurnal  variation.  All  data  collected  were  used  for  boundary  specifi¬ 
cation,  but  diurnal  variations  were  not  important  and  need  not  have  been 
collected. 

90.  The  13-hr  surveys  conducted  at  the  downstream  boundary  provided 
adequate  data  for  characterizing  intra-  and  intertidal  variation  of  tempera¬ 
ture,  salinity,  and  DO.  Examination  of  time  plots  of  BOD,  NH^,  and  P  revealed 
that  the  boundary  constituent  values  could  have  been  averaged  for  each  date 
and  linearly  interpolated  between  dates. 

Initial  conditions 

91.  Extended  simulation  of  salinity  starting  with  a  concentration  of 
0.0  ppt  revealed  that  equilibration  was  achieved  in  10  to  20  tidal  cycles. 

The  exercise  revealed  that  "hot  starts"  following  simulations  extending  for 
greater  than  20  tidal  cycles  were  adequate  for  the  specification  of  initial 
conditions . 


Model  adjustment 

92.  Errors  in  model  geometry  manifested  themselves  as  tidal  phase  and 
water  velocity  errors.  It  is  hypothesized  that  inclusion  of  overbank  flooding 
would  substantially  eliminate  the  phase  and  tidal  errors;  adjustment  of  the 
Chezy  coefficient  had  little  effect. 


93.  Discrepancies  between  simulated  and  measured  salinity  may  be  due  to 
errors  in  advective  transport.  The  simulated  surface  velocities  were  less 
than  measured,  and  adjustment  of  the  horizontal  dispersion  coefficient  had  no 
effect  on  simulated  salinity  concentrations. 

94.  Simulated  and  measured  profiles  of  temperature,  BOD,  and  DO  were 
nearly  equivalent.  The  results  were  quite  good  considering  the  difficulty  in 
specifying  the  head  boundary  conditions  and  matching  the  advective  transport. 
Experimental  results 

95.  Experiments  that  were  conducted  demonstrated  the  utility  of 
CE-QUAL-VJ2  for  investigating  changes  in  upstream  regulation  and  evaluating 
permit  applications. 

96.  Two  sensitivity  experiments  were  conducted  in  order  to  evaluate  the 
effects  of  imprecise  specification  of  head  boundary  concentrations.  The  first 
consisted  of  a  factorial  arrangement  of  DO  and  BOD  head  boundary  concentra¬ 
tions  using  simulated  tidally  averaged  DO  concentration  as  the  response 
variable.  The  results  indicated  that  increasing  the  head  boundary  BOD  con¬ 
centration  100  percent  decreased  the  tidally  averaged  DO  concentration  1  per¬ 
cent;  decreasing  the  head  boundary  DO  concentration  50  percent  decreased  the 
tidally  averaged  DO  concentration  10  percent. 

97.  The  second  sensitivity  experiment  compared  the  relative  Importance 
of  the  head  boundary  algal  concentration  with  the  algal  maximum  growth  rate. 

A  100  percent  error  in  estimation  of  algal  head  concentration  would  result  in 
an  average  algal  mass  error  of  65  percent,  whereas  a  100  percent  error  in 
estimation  of  maximum  algal  growth  rate  would  result  in  an  average  algal  mass 
error  of  11  percent. 


Recommendat ions 


98.  Detailed  cross-sectional  profiles  of  the  study  area  need  to  be 
obtained.  The  cross-sectional  profiles  should  Include  areas  that  are  tidally 
inundated . 

99.  Specification  of  external  head  boundary  profiles  was  critical.  The 
external  head  boundary  profiles  may  be  more  easily  specified  at  the  mouth  of 
the  estuary,  where  vertical  homogeneity  can  frequently  be  assumed. 

100.  The  upestuary  flow  boundary  should  be  specified  above  the  fall 
line  in  order  to  avoid  tidal  Interaction. 
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101.  As  a  starting  point,  it  is  recommended  that  the  following  water 
quality  constituents  be  sampled: 

a.  Inorganic  suspended  solids. 

b.  Salinity. 

c.  Dissolved  organic  matter. 

d.  Particulate,  nonliving,  organic  matter. 

e.  Algae. 

£.  Inorganic  phosphorus. 

£.  Ammonia  nitrogen. 

h.  Nitrite  plus  nitrate  nitrogen. 

i.  Temperature. 

Dissolved  oxygen. 

k.  Sediment  oxygen  demand. 
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